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ABSTRACT 

Context. The combination of optical and near-infrared (NIR) colours has the potential to break the age/met allicity 
degeneracy and offers a better metallicity sensitivity than optical colours alone. Previous tudies of extragalactic globular 
clusters (GCs) with this colour combination, however, have suffered from small samples or have been restricted to a 
few galaxies. 

Aims. We compile a homogeneous and representative sample of GC systems with multi-band photometry to be used in 
subsequent papers where ages and metallicity distributions will be studied. 

Methods. We acquired deep K-band images of 14 bright nearby early-type galaxies. The images were obtained with the 
LIRIS near-infrared spectrograph and imager at the William Herschel Telescope (WHT) and combined with optical 
ACS g and z images from the Hubble Space Telescope public archive. 

Results. For the first time GC photometry of 14 galaxies are observed and reduced homogeneously in this wavelength 
regime. We achieved a limiting magnitude of K ~ 20 — 21. For the majority of the galaxies we detect about 70 GCs 
each. NGC4486 and NGC4649, the cluster-richest galaxies in the sample contain 301 and 167 GCs, respectively. We 
present tables containing coordinates, photometry and sizes of the GCs available. 
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1. Introduction 

Understanding how galaxies formed and evolved is one 
of the fundamental problems in modern astrophysics. 
Spheroidal systems (ellipicals, lenticulars and spiral bulges) 
contain ~ 75 % of the Local Universe's stellar content 
(Re nzini||2006[ ), and much of our knowledge on galaxy for- 
mation and evolution comes from detailed studies of galax- 
ies in this close portion of the Universe. In the hierarchi- 
cal merging framework of structure formation, early-type 
galaxies are the products of mergers, and star formation 
is governe d by the amount of ga s involved in the merger 
event (eg. Kauffman et al. 1993). This model is very suc- 
cessful in explaining structure for mation, being the prod - 
uct of the ACDM cosmology (eg. |Komatsu et a l. 2009), 
which is the best at describing the observations of the cos- 
mic microwave background, large scale structure and the 
accelerated expansion of the Universe with the Big Bang 
theory. Nonetheless, it remains to be understood how bary- 
onic galaxy formation occurs within this framework. 

Globular clusters (GCs) are found around all major 
galaxies, and young massive clusters (YMCs) ar e being ob- 
served in many mergers and standard spirals (Whitmore 
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* Table 4 is only available in electronic form at the CDS 



via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 



http://cdsweb.u-strasbg.tr/cgi-bin/qcat'.'' J/ A+A/ 



& Schweizer 
alaxies (ear 



1995) [Larsen fc Richtler| [l999). If (massive) 
y-type) formed through mergers, these YMCs 



are good candidates to be the progenitors of the GCs we 
observe today in many early-type galaxies. This ubiquity 
of star clusters around all major galaxies suggests that GC 
formation traces important events in the history of their 
parent galaxies. Even though we cannot state anymore that 
all GCs are simple stellar populations (SSPs) due to t he evi- 
dence of multiple populations in several Galactic (eg. |Bedin| 
et al.||2004| |Piotto et al.||2007[) and Large Magell anic cloud 
el 



star clusters ^antiago et al.||2002| [Mackey|2008[ ), they 
still more easily modelled than the mix of stellar pop 



are 
popula- 



tions that make up the integrated light of a galaxy. 

A thorough understanding of the evolutionary histories 
of galaxies must involve both their GC systems and the 
integrated light. Integrated light studies suffer from a de- 
generacy between age, metallicity and burst strength, which 
is difficult to overcome. GCs on the other hand, have the 
advantage of being discrete and separable star formation 
tracers. Generally, a significant population of GCs with a 
particular age would be a strong indication that the host 
galaxy experienced a burst of star formation at the corre- 
sponding epoch. If such a burst occurred recently, within 
the last few Gyrs, and involved a significant fraction of the 
stars, it should also be detectable in the integrated light. A 
significant fraction, if not the vast majorit y of all the stars 
are formed in clusters ( Lada & Lada 2003 ) and a consider- 
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able part of these do not remain bound for a Hubble time 
and probably help build up the integrated light of a galaxy. 
Therefore, a galaxy's integrated light should be correlated 
in many ways to that of its GC system. Studies of GCs and 
their host galaxies are thus complementary and evidence of 
star formation seen in the integrated light should also leave 
imprints in the GC system. 

The colour distribution of GC systems in most galaxies 
is bi/multimodal in optical colours. Many stud ies, espe- 



number statistics, metallicity distributions of GC systems 
in early typ e galaxies are poo rly constrained spectroscopi- 
cally. M 87 ( |Cohen et al.|1998[ ) shows marginal evidence for 
a bimodal as opposed to unimodal metallicity distribution. 
Nevertheless, there is evidence that this could also be due 



to the blue tilt or lumino sit y-colour relation (see eg. Strader 



et al. |2006l |Harris|2006l) a^ 



ciall y spectroscopic ones (eg. Strader et al- 
and |Strader et al.|2005|[Cenarro et al.|2007 for indirect ev- 
idence), show that this feature is due to 2 old populations 
(> 10 Gyr) that differ in metallicity. Bimodal colour (metal- 
licity) distributions may indicate distinct cluster formation 
mechanisms for the two populations, with age differences 
of a few Gyrs (~ 2 Gyrs) allowed within the unc ertainties. 
Moreover the re is sufficient evidence available (Brodie & 
Strader|2006 ) to place metal-rich GC formation along with 
the bulk of the field stars in their parent galaxies. Most GC 
formation scenarios introduced to account for co lour bi- 
modality (lAshman & Zepf|1992l|Forbes et al.|1997 



Strader et al. ( 2007 ) reanalysed 
thelCotoiet^llgOOSI) sam ple of GCsm NUG4 472, which 
does not show the blue tilt ( Strader et al.||200"6 ) and found 



2007 for direct bimodal metallicity distributions. 



al.|1998[|Beasley et al.|2002[|Strader et al.|200T^ and |Rhod~ 



Cote et 



et al.||2005[ ) assume different formation/assembly channels 
for the red and blue clusters although these clusters have 
very similar interna l prop erties such as masses and sizes. 



Muratov & Gnedin 
for the formation of 



2010| propose a common mechanism 
oth sub-populations in which the bi- 
modal metallicity distribution arises naturally in many of 
the model realisations from the combination of the history 
of galaxy assembly and the amount of cold gas in proto- 
galactic systems. This scenario is built upon the hierarchical 
merging framework, in which early mergers of smaller hosts 
create only blue clusters, and subsequent mergers o f pro- 
genitor gala xies produce both blue and red clusters. |Yoon| 



et al. ( 2006[ ), however, challenged the metallicity bimodal- 
ity interpretation. These authors claim that the horizontal 
branch stars (HB) bring up a non-linear behavior on the 
metallicity-colour relation of GCs. This non-linearity is in 
turn responsible for making a coeval group of clusters with 
a unimodal metallicity spread exhibit colour bimodality. 

We want to understand metallicity distributions of GC 
systems and how their age distributions are shaped. The 
most direct way to measure the metallicity and derive ages 
for GC systems would be to take spectra of its individual 
members. Spectroscopy for a large number of GCs is very 
time consuming; it is difficult to collect data for represen- 



going effort of the extragalactic GC community (Cohen et 


al. 


1998 


Strader et al. 2005, Puzia et al.|2005b[|Cenarro et 


al. 


2007 


) to collect a significant amount of data for a large 



number of galaxies. The majori ty of spectroscopic s tudies 



reach around 2 dozen clusters (Strader et al. 2005 Puzia 
et al.|[2005b[ ) biased to the brightest parts. It still remains 
unclear if these are representativ e of the system as a whole, 
with the exception of a few cases ( Puzia et al. 2005a on M 31 
reached ~70 GCs ; |Woodley et al.||2009| on NGU5128 and 
Cohen et al.||1998 on M87 who reached over 100 clusters). 

The general result is that the studied clusters are pre- 
dominantly old (> 10 Gyrs) but there is some disagreement 
on the fracti on of GCs with ages of 5 — 6 Gyrs. While some 



studies (eg. Strader et al. 2005 Cohen et al. 1998) find 



little ev idence for young GCs, t he study by |Puzia et al 



(2005b) finds that up to 1/3 of their 143 clusters in differ- 
ent E/S0 galaxies are younger than 10 Gyrs. Due to small 



Given the difficulty of current spectroscopic studies in 
reaching representative statistics, NIR imaging appears as 
a powerful alternative. One can achieve much better statis- 
tics, reaching fainter GCs than the deepest spectroscopic 
studies with shorter integration times. It has been shown 
that the colours that best represent the true metallicity 
distributions are the combination of optical an d NIR (eg. 
iPuzia et al. 



2002 

NIR depends mam 



Cantiello fc Blakeslee||2007[ ). While the 
y on the red giant branch being metallic- 



ity sensitive, the optical depends on stars near the main se- 
quence turn-off point being both metallicity and age sensi- 
tive. By using the optical/NIR combination we have, there- 
fore, a better chance to break the age-metallicity degener- 
acy inherent in optical colours alone. Hence, studying GCs 
in the NIR is an alternative to spectroscopy in revealing 
their actual metallicity distributions. 

There have been several attempts to study GC systems 



with NIR imaging with a variety of instruments (Hempel 
et al.||2007a[ |Hempel et al.||2007b[ |Kotulla et aLl|2008| and 



Brodie fc Strader|2006 for a compilation). With the excep 
tion of a few cases, the detectors had generally small fields 
of view (~ 2 x 2 arcmin 2 ). In total around a dozen sys- 
tems have been studied with the optical/NIR technique. 
The general result is that most of the studied galaxies 
host predominantly old clusters, but there is still a pos- 
sible significant fraction of intermediate age GCs in some, 
with the strongest and most extensively studied cases be- 
ing NGC1316, NGC4365 and NGC5846. While the first 
one has a power-law rather than lo g-normal GC luminosit y 



function for the metal rich GCs (Goudfrooij et al. 2001) 



which is strong evidence for a younger population, the GC 
systems of the other two do not present any outstanding 



Yamada et al. 


2006) 


the optical/ N 


IR ap 



proach there has been only one case so far. |Kundu &: Zepf| 



( |2007[ ) studied the (I-H) colour distribution for 80 GCs in 
M87 with NICMOS and WFPC2 data and found a clear 
bimodal distribution. 

Improvements in NIR SSP models over th e p ast few 
years have been significant (eg. Maraston|2005 and Marigo 



|et al.| [20081 for the TP-AGB perclval et ah] |2009| for 
a— enhancement estimates). Along with this progress, the 
full potential of the optical/NIR technique should be fur- 
ther explored. The purpose of this project is to establish 
a large and homogeneous optical/near- infrared data set of 
GC systems in elliptical and lenticular galaxies where dif- 
ferential comparisons among the different GC systems are 
possible. In the present paper the survey is introduced, the 
data reduction procedures are thoroughly explained and 
catalogues of the GC systems with g, z and K photometry 
and sizes are made available. Ages and metallicity distribu- 
tions will be studied in subsequent papers. 
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2. Galaxy sample selection 

The target galaxies were se lected from the SAURON sam- 
ple ( |de Zeeuw et al. 2002). Of the 72 galaxies that make 
up the SAURON survey, we added three more cases that 
were observed with the instrument and restricted our- 
selves to 16 of these 75 galaxies that were observed with 
HST/ACS in g (F475W) and z (F850LP) to keep the sam- 
ple homogeneous. Of these 16, we obtained K-band imag- 
ing for 14. Two of the 16 were not observed due to bad 
weather conditions and override programs during the ob- 
serving runs. The SAURON data provide a wide range of 
kinematic and stellar population parameters for the galax- 
ies, including line indices sensitive to age and metallicity 
(H/3, Fe5015, Mg6) and their gradients, presence of dis- 
tinct kinematic components and rotation parameters. We 
required distance moduli (to — M) < 32 and absolute mag- 
nitudes Mb < —19 to ensure the detection of significant 
numbers of GCs. While this sample is not representative 
of early-type galaxies as a whole because we do not cover 
low-mass galaxies (Mb > —19), it is quite representative 
within the SAURON sample, including galaxies with either 
boxy or disky isophotes, being slow or fast rotators, and 
having a variety of kinematic properties. The characteris- 
tics of these galaxies are listed in Tables [l] and [2] where we 
also list esti mates of ages and metallicities from integrated 
light studies flSanchez-Blazquez et al.||2006| |Yamada et al. 
2006 Kuntschncr et al.||2010D and the estimated number 



of GCsand their specificTrequ ency (Sat) ( Peng et aL]|2008 
and |Kundu fe Whitmore|[200l| ) . 

Our sample contains 10 galaxies classified as ellipticals 
and 4 as SO's; 5 are slow rotators and 9 are fast rotators. All 
but two of the galaxies (NGC 3377 and NGC 4278) belong 
to the Virgo Cluster of Galaxies. In the following a short 
description of earlier work on the 14 GC systems is pre- 
sented, and the most important results on kinematics and 
integrated stellar populations from the SAURON survey 
are summarised. 

NGC 3377: The GC system of thi s E5/E6 type galaxy, 
locate d in Leo Group I, was studied by |Kundu fc W hitmorc 
(2001) with (V-I), and a likely bimodal colour distrib ution 
was found. It was suggested by Kormendy et al. ( |1998[ ) that 
this galaxy contains a black hole. Contains spatially ex- 
tendecflsig ns of SSP equivalent^] younger ages (Kuntschner 
et~al|2 010). 

NGC 4278: |Forbes| ( fl996| studied the GC luminosity 
function of this E1/ E2 galaxy and confir med its location 
in the Coma I cloud. Forbes et al. ( 1996 ) studied the (V- 
I) colour distribution of th e GC system finding no ob- 
vious bimodality. Later on, Kundu & Whitmore (2001) 
with (V-I) , found a likely bim odal colour distribution. Thi s 
galaxy hosts a LINER AGN ( |Veron-Cetty fc Veron|2006 1 . 
Kuntschner et al72 010 report ages consistent with old stel- 
lar populations. 

The following galaxies are also part of the ACS Virgo 
Clust er Survey (ACSVCS) ( [Jordan et al1|2004l |Peng et al. 
2006 ) , which made use of the (g-z) colour. Unless otherwise 
stated, this is the colour used in the GC systems studies. We 

1 spatially extended, i.e covering la rge parts of the observed 
field of view dKuntschner et al.|2010h. 

9 . ' 1 ! ! A 



2 SSP equivalent ages, ie., when any observable stellar popu- 
lation parameter is the one the object would hav e if it formed 



at a s ingle age with a single chemical composition ( Trager et al 
2008) 



refer the reader to the appendix of Ferrarese et al. ( 2006 ) 
for notes on the morphology of the individual galaxies. 
NG C 4365: The GC sy stem of this E3 galaxy was stud- 



Whitmore 
vious bimo 



ied by Forbes et al. ( 1996) and subseque ntly by Kundu & 
( 200 ip and |Larsen et al.| (|2001|) who fou nd no ob 
idality in the (V-I 



ity in the (V-I) colour distribution. Peng et 
al. (2006) reported a bimodal distribution with a broad red 
peak and an excess of red GCs when compared to similar lu 
minosity ACSVCS galaxies. This galax y has also been stud - 
ied with optical / NIR 



(2005) 



photometry by Puzia et al.| (|2002 



by |La rsen et al. (2005) with V, I and K, by Ku ndu et a. 
with g, I and E 



al. 



and spectr oscopically by Larsen et 



|al.| (|2003) and Brodie et al. ( 2005[ ). Some of these authors 
find evidence for a significant fraction of young GCs. It 



hosts a kinematically distinct component (KDC) (Bender 
1 998[) and has ages consistent with old stellar populations 
( Davies et aL|2001[ ) 



NGC 4374: Also known as M 84, the GC sy stem of this 
galaxy was studied by Gomez fc Richtler ( 2004 ) with (B-R) 
and by |Peng et al. ( 2006 ) , who find a bimodal colour dis- 
tribution. This last study finds a sparsely populated, blue 
peak when compared to same luminosity- t ype ACSVCS 
early- type galaxies. It is a Seyfer t 2 AG N ( |Veron-Cetty 



& Veron 2006). Kuntschner et al. (20101 report ages con 



sistent with old stellar populations 



NGC 4382: |Peng eTal] ( p006| ) report the GC system 
of this galaxy, also known as M85, to be the best case for 
a trimodal colour distribution among the ACSVCS-studied 
galaxies. It contains a spatially constrained, central (po st) 
starbursl^] ( |Kuntschner et al.||2010| |Shapiro et"aT1|2010| ). 

NGC 4406: This S0/E3 galaxy, also known as M86, 
is a noteworthy gala xy with on-going dust stripping 



( Ferrarese et al. 
Forbes 



et al 



2006J). It is also part of the studies by 
}l$9ty and jLarsen et al.| (|2001|), with no 



strong bimodality being reported. Both studies w ere car- 



ried out with (V-I) . Similar results were found by Kundu 
fc Whitmore (2001[) (with V-I photometry) and Peng et al. 



(200 6THT was also studied with B, V and R photometry. 



by Rhode fc Zepf] ( |2004[ ) who reported a mu ltimodal colour 

distribution. It hostsT KDC (|Bender||1998|). 

NGC 4473: While |Kundu fc Whitmorel (|200ll (with 



V-I photometry) and Peng et al. ( 2006D report a b i modal 



colour distribution for this system, Larsen et al. (2001) 



(wi th V-I photometry) find marginally significant bimodal 



ity. |Kuntschner et al. (20101 report ages consistent with old 
stellar populations. 

NGC 4486: The Virgo central giant elliptical, also 
know as M87, is the most extensively studied galaxy in the 
sample, hosting the largest numbers of clu s ters (^15000, 



Peng et al ||2008| ). Elson fc Santiago| (|1996[), |Larsen et al 



( 200ip both with (V-I) and |Peng et al.| ( |2006l ) reporTsig 



nificant colour bimodality. The GC system of this galaxy 
has been reported to host a correlation between colour and 
luminosity for indiv idual, metal-poor GCs, the blue tilt 



( Str ader et al 
I Peng et al 



tL"1|2006| both with (g-z), |Mieske et al 
2009[ |Harris||2lX)9| with g',r' an d i). Ev 



2006 
Evidence 



for signific ant, intermediate- age population ( Kissler-Patig 
et al.|200"2 ) is not found through optical/NIR i maging with 



V, 1 and K photometry. Kundu & Zepf ( 2007 ) find a clear 
optical/NIR (I-H) bimodal distribution. Cohen et al.|(|1998 ) 



spatially constrained central starburst, i.e, the star forma- 



tion is connected to a thin embedded disk (Kuntschner et al. 
20101. 
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present a large spectroscopic sample of GCs in M87, and 
they are inferred to be uniformly old (> 10 Gyrs). A clear 
bimodality in [Fe /H] and Mgb is not seen in that study. It 
is a L I NER AGN ( |Veron-Cetty fc Veron|2006| ). |Kuntschner 
et al. (2010) report ages consistent with old stellar popula- 



tions. 



NGC 4526: |Peng et al.| ( |2006| ) report significant colour 
bimodality. Contains residual, ongoing star formation in 



a central disk/ring structure (Kuntschner et al. 20101 



Shapiro et al. (2010) estimate the star formation rate for 

l 



this galaxy to be 0.37 Mq yr 

NGC 4552: The GC system of this galaxy, also k nown 
as M89, was s tudied by |Kundu fc Whitmore| fl200lD and 
Larsen et al. (2001) with (V-lj colour and by Feng et 
aT ( |2006[ ) who report likely colour bimod ality. It is a 
Seyfer t 2 AGN (|Veron-Cetty fc Veron]|2006D and hosts a 
KDC ( Emsellem et al.|2007p . A pronounced Hj3 absorption 
turn-up is seen in this galaxy, although the presence of a 
spatially constrained, ce ntral (post) starburst is uncertain 
( Kuntschner et al |2010j ). 

NGC 4570: Like NGC 4374 the colour d istribution of 
this galaxy is likely bimodal ( Peng et al.|2006 ), with a small 
blue peak if compar ed to similar luminos ity-type ACSVCS 
early-type galaxies. Kotulla et al. (2008) find a significant 
population of intermediate-age g lobular clusters wi t h g, z 
and K optical/NIR photometry. Kuntschner et ah] (2010) 
report ages consistent with old stellar populations 



NGC 4 621: | Kundu fc Whitmore| pOlj ) with (V-I) 
colour and Feng et al. P006|) report obvious bim odal- 
ity. This galaxy hosts a KDC ( |Emsellem et al. 2007) and 
Kuntschner et al. (2010) report ages consistent with old 
stellar populations. 



NGC 4649 : |Kundu fc Whitmore| ( 2001|) an d |Larsen et 
ah] ( |2001[ ) with (V-l) and |Peng et al.| pioop r eport sig- 
nmcant colour bimodality. Forbes et al.| ( |2004p find that 
the underlying galaxy starlight has a similar density profile 
slope and (g-i) colour to the red GCs. The GC system of 
this galaxy, al so know as M60, has been reported to exhibit 



the blue tilt (Strader et al. 2006 Mieske et al. 2006 Lee 



et al. 


2008|). Some GCs have also been observed spectro- 


scopically (Bridges et al. 2006 


Pierce et al. 20061 


Hwang et 


al. 2008 


|. Bridges et al. (|2006 


concluded tnat th 


e GC sys- 



from its kinematics. Pierce et al. ( |2006[ ) find most GCs to 
be old (> 10 Gyr), but several intermediate-ages are present 
(2 — 3 Gyrs). These authors also find a trend of decreasing 



a-elem ent ratio with increasing metallicity. |Hwang et al. 
( |2008[ ), on the other hand, report that the NGC4649 GC 
system exhibits significant overall rotation and also a dar k- 



2003) 



matter halo. It hosts a black hole (Gebhardt et al 

NGC 4660: The GC sy stem of this galaxy l^oun cTto 
be not significan tly bimodal (|Kundu fc Whitmore 2001|with 
the (V-I) colour, |Peng et al.|2006p . |Kuntschner et al.|Q2010| 
report ages consistent with old stellar populations. 

It appears that some studies find colour bimodal- 
ity while others do not, even for the same galaxies (eg. 
NGC 4278). These differences can be attributed to the dif- 
ferent colours used in the different studies, which sample 
stars in different evolutionary stages. Another factor for 
different shapes in the colour distributions are the different 
data sets used in different studies. Different data sets have 
different photometric uncertainties with large uncertainties 
likely causing a bimodal distribution to become unimodal 
if the range in wavelength is too narrow. 
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3. Observations and data reduction 

3.1. LIRIS K-band imaging 

Imaging in the Ks-band (from now on referred to as K) 
was obtained during 2007, 2008 and 200 9 with the LIRIS 
ne ar- infrared spectrograph and imager ( |Acosta Pulido et 
aT1 |2003[ panchado et al ||2004| ) on the William Herschel 
Telescope (WHTJ in La Palma, Spain. A summary of the 
LIRIS observations is given in Table [3] This instrument has 
a field of view of 4.2' X 4.2' and a plate scale of 0.25" pix^ 1 . 
We adopted a dithering pattern that consisted of taking ex- 
posures in 5 different positions across the detector, starting 
at the centre, moving to the upper right quadrant and in 
a clock wise rotation pattern spanning all the 4 quadrants 
with a jitter from 12 " to 20 ". We took 4 exposures in each 
of the 5 positions of the detector. Each exposure consisted of 
15 seconds to cope with the high count rate from the back- 
ground in the K-band. The total exposure time reached for 
each galaxy is also given in Table [3j 

Even though the fields are dominated by the host 
galaxy, we did not adopt the traditional method of sky sub- 
traction in the NIR for crowded fields. In that method one 
takes a sky exposure for every science exposure and spends 
the same amount of time observing the sky as observing 
the actual target. We only took a few sky exposures (~ 25 
frames per galaxy), and these few frames allowed a fairly 
good sky subtraction on a subset of the sky exposures and 
a good enough ellipse modelling and removal of the galaxy 
light. The images resulting from this operation were then 
used to create the actual "sky" image. 

For the 2007 run, 8 standard stars were observed during 
the night of 27/02 and 7 during the night of 28/02, and both 
nights were found to be photometric. For the 2008 run, 7 
standard stars were observed during the night of 18/03, the 
only photometric night. For the 2009 run, 6 standard stars 
were observed during the night of 15/03, the only night 
considered to be photometric. During the last run in 2009, 
we had problems with Calima dust storms from the Sahara 
desert, which severely affected the data taken on the night 
of March 13. 



3.2. LIRIS data reduction and calibration 

The data reduction was performed with the standard pack- 
ages in IRAF and the LIRISDR package written by Jose 
Acosta Pulido and kindly provided to us by the author. The 
shifts between exposures were determined by measuring the 
centroids of 2 (min) to 4 (when possible) isolated bright 
stars in each exposure. The first step, was the correction of 
the pixel-mapping anomaly of LIRIS with LCPIXMAP in 
LIRISDR before any reduction procedure took place. When 
the two-dimensional image is reconstructed from a one di- 
mensional array, some of the pixels in the LIRIS image end 
up in places where they should not be. This affects mainly 
the lower left of the first quadrant ( Schirmerj|2007 ). 

A crude sky subtraction using sky frames within 
< 30 min from the science exposures was performed. The 
resulting images from this crude sky subtraction were then 
combined after flat-field correction in each individual ex- 
posure. From this combined image an ellipse model was 
created with the STSDAS tasks ELLIPSE and BMODEL. 
This ellipse model was then multiplied by the flat field, 
shifted and subtracted from each original exposure. We 



were left with exposures without the central galaxy cusp 
on them. These " galaxy- free" images were median com- 
bined following the standard procedure of sky subtraction 
in the near- infrared: a "running mean" of the 4 exposures 
obtained closest in time (2 before and 2 after) was con- 
structed for each frame. The result of this operation yielded 
the sky image to be subtracted from each original exposure. 
This technique of sky subtraction has optimal results for 
the simplest E0 galaxies, but as the galaxy morphology be- 
comes more complex, the ellipse model fit becomes poorer. 
Fairly good ellipse models are obtained for all types of E- 
type galaxies but poor models are achieved for SO galax- 
ies as these galaxies contain an additional disc component. 
Since ELLIPSE is not able to model such a component, in 
the disc-dominated part of the galaxy our photometry will 
be compromised. 

In Fig. [T] we illustrate some of the steps of the sky sub- 
traction. Panel (a) shows a 15-second science frame that 
is sky subtracted by the few science frames (taken within 
< 30 min from it) and results in panel (b). A stack of im- 
ages like panel (b) are averaged together (60 — 200 frames), 
and the resulting image is used to create an ellipse model 
fit, shown in panel (c). If one subtracts an image like (c) 
from an image like (a), the result is (d). Median combining 
several images like (d) with different dither positions will 
allow for the removal of its model residual, seen in the cen- 
tre of this panel as well as any bright foreground star or 
background galaxy that is in the image. The result of this 
combination will be an image like the one shown in panel 
(e), an average sky frame. Subtracting the constructed sky 
from the original science frame (a) will result in a sky sub- 
tracted science frame, shown in panel (f). 

After flat-fielding the sky-subtracted images we cor- 
rected for the vertical gradient with LICVGRAD in 
LIRISDR, another anomaly suffered by LIRIS. Some ex- 
posures present a discontinuous jump between the upper 
and the lower two read-out quadrants, with the amplitude 
of this jump depending o n the order of a n image in a series 
of continuous exposures ( Schirmer|2007 1 . This happens be- 
cause when the telescope is moved from one dither position 
to another it takes about 3 exposures for the detector to 
reach equilibrium. Each image is then corrected individu- 
ally. Following the vertical gradient correction, the images 
were shifted and averaged resulting in a reduced, final im- 
age. Figure [2] shows these images for our galaxy sample. 
Another elliptical model was created for each galaxy and 
subtracted from the final, reduced image. We did not use 
the same model constructed through a crude sky subtrac- 
tion (eg. Fig. [I] panel (c)) as it is based on a small fraction 
of the images, and the final, reduced images contain ~ 800 
exposures. In this model-subtracted image we performed 
photometry with PHOT after object detection in the opti- 
cal ACS images, as described in Sect. 3.3. 

To calibrate our data photometrically, we assumed a 
transformation of the following form: 



+ z k -p k x x ka 



(1) 



where K s and k s are the standard and instrumental mag- 
nitudes, and pk and the atmospheric extinction 
coefficient and the effective airmass, respectively. To deter- 
mine pk we needed a broad range in airmass values, and the 
standard star exposures didn't span the needed range. We 
then used both standard stars and individual galaxy expo- 
sures that contained bright stars taken at different times 
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Fig. 1. Ilustration of some of the important steps involved in the sky subtraction for NGC4374. (a): science frame, (b): 
rough sky subtracted science frame, (c): ellipse model, (d): ellipse model subtracted frame, (e): an average sky frame, 
(f): sky subtracted science frame 



during the night, ie., covering sufficient airmass values to 
derive p k during the 2007 run. An estimate of p k was calcu- 
lated to be the slope of the linear relation between airmass 
and the observed magnitudes of the same standard star 
and/or bright stars in the individual science frame expo- 
sures taken during the night. We adopted a value of pk = 
0.07 ± 0.02 (unit airmass) -1 , determined as the mean of the 
slopes found for different bright and/or standard stars in 
these exposures. There were no previously documented val- 
ues for pk at this wavelength for La Palma. Typical values 
found in the literature for p k for sites such as C erro Paranal 
and La Silla are p k = -05 (unit airmass) -1 (Puzia et al. 

0.08 (unit airmass) -1 for 



2002 



Larsen et a l. 2005) and p k 



Las Campanas (Hempel et al. 2007a). To determine x kcii , 
the mean values for combined exposures, we used the fol- 
lowing expression, which accounts for the airmass in each 
of the individual exposures and where the subscript i des- 
ignates the i th exposure 



Xcti 



0.4 x p k 



(2) 



Equation (2) is found by summing up the individual con- 
tributions of each airmass value in the equality, and solving 
for x c ff 



Fn 



F k x 10 



-0.4p fc xa;eff 



(3) 



Following the standard star observations, we determined 
the zero point using an aperture of 20 pixels. For the 2007 
run it was found to be z k = 23.450± 0.004, for the 2008 run 
z k = 23.110±0.004 and for the 2009 run z k = 23.070±0.008. 
For the 2008 and 2009 runs only the first night was found 
to be photometric and an extra step, ie. the photometric 
correction, was added to the equation. These values are 
listed in Table [3j The errors in the z k values are the stan- 
dard errors of the mean. Objects were detected in the ACS 
images (see Sect. 3.3) and transformed to the LIRIS coor- 
dinate system. Aperture photometry was performed with 
PHOT in a 5 pixel aperture in the final, model-subtracted 
image. We estimated photometric corrections by measur- 
ing the magnitudes of some bright objects (~ 10) for each 
galaxy in a combined image with all exposures and only 
with the photometric night exposures. The weighted mean 
difference between these for each galaxy was taken to be the 
photometric correction adopted value. An effort was made 
to always obtain some exposures of each galaxy under pho- 
tometric conditions. The only galaxy where this was not 
possible was NGC 4278 with exposures taken on the nights 
of March 22/03/08 and 24/03/08. An analysis of the rela- 
tion between airmass and instrumental magnitude for a few 
bright objetcs in all exposures showed that even though the 
latter night was not photometric, its exposures were not as 
severely affect by extinction as for the first night. An es- 
timate for the extinction, through exposures taken in the 
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Fig. 2. Liris final reduced K-band images for the 14 galaxies. North is up and east to the left. Each image is ~ 270" on 



a side, corresponding to ~ 24Kpc for Virgo Cluster galaxies, 
colour bar is given in units of counts. 



21 Kpc for NGC 4278 and - 15 Kpc for NGC 3377. The 



first night, was determined to be 0.030 ± 0.006 mag. Since 
this value is small and we are not certain about the extinc- 
tion for the last night, we did not apply the photometric 
correction for this galaxy. As can be seen in Table [3j the 
uncertainties in the photometric corrections are < 0.007 
for the second run and < 0.03 for the third one. Aperture 
corrections were calculated with APPFILE from 5 to 20 
pixels and cross-checked by determining the difference be- 



tween the magnitudes measured with a 5 pixel aperture to 
a 20 pixel one, and the values proved to be consistent. With 
the uncertainties on the zero points, photometric and aper- 
ture corrections listed above, we estimate the calibration to 
be accurate to - 0.04 for M87, ~ 0.03 for NGC 3377 and 
~ 0.05 for NGC 4660. We show these estimates as examples 
for each of the observation runs. These accuracy estimates 
were calculated by adding in quadrature the uncertainties 



10 



Chies-Santos et al.: GCS of early- type galaxies in the NIR 



in the zero points, photometric and aperture corrections. 
The magnitudes were correct ed for galactic extinc tion us- 
ing the values in NEER from |Schlegel et~aT| < |1998| >. These 
quantities are also listed in Table |3| 



3.3. H ST ACS imaging 

ACS pipeline processed images of the 14 galaxies were 
downloaded from the MAST archive. The data for 12 of 
them were obtain ed as part of the ACS Virgo Cluster 
Survey (ACSVCS; |C6te et al |[2004| ) and use the F475W 
(~Sloan g) and F850LP (~Sloan z) filters. The exposure 
times were 750 s and 1210 s in F475W and F850LP, respec- 
tively, split into 2 and 3 exposures for cosmic-ray rejection. 
From now on we refer to these filters as g and z. For the 
galaxies not located in the Virgo Galaxy Cluster, NGC 3377 
and NGC 4278, there were 1 and 2 pointings, respectively, 
that covered the LIRIS field of view. The ACS Program 
IDs are 10435 and 10835. For the former the exposure times 
were 1380 s for g and 3005 s for z, divided into 4 exposures 
for both g and z. For the latter galaxy, the exposure time 
for the z band was 1200s for the two pointings, split into 3 
exposures; for the g band the central pointing has 676.2 s 
of exposure time and for a west pointing 702 s split into 2 
exposures. 

The data were obtained with the wide field channel 
on ACS, which has a pixel scale of 0.050" pix -1 and a 
field of view of 3.36' x 3.36'. This is roughly the field 
of view of the final, combined LIRIS images considering 
the dithering pattern. This gives us a radial coverage of 
- 8 kpc for NGC 3377, ~ 11 kpc for NGC 4278 and ~ 12 kpc 
for the Virgo Cluster galaxies. In terms of the effective 
radius (R e ) of the galaxies, it is 0.29i? e for NGC 3377, 
0.22 R e for NGC 4278 and for the Virgo Cluster galax- 
ies 0.08 R e for NGC 4486 and 1.3 Re for NGC 4660. They 
were then combined with the MULTIDRIZZLE task in the 
STSDAS. DRIZZLE package, which removes the geometric 
distortion and normalises the images to exposure times of 
1 s. A 20 x 20 median filter image was subtracted from the 
multidrizzled image and the sources were detected on the 
resulting image with DAOPHOT in the z band image. Since 
there is significantly more background noise near the cen- 
tre of galaxies than farther out, more spurious objects are 
detected in the central part of the galaxy. At this point 
nothing was done to remove these spurious objects as they 
are expected to be automatically removed from the sam- 
ple with the application of the criteria to refine the GC 
selection, explained in Sect. 4. Aperture photometry was 
performed in a 5-pixel aperture. The following aperture cor- 
rections were used for the Virgo Cluster galaxies: —0.09 in g 
and —0.15 in z from 5 to 10 pixels. These values were taken 
from Strader et al. (2006), who derived them as median 
corrections from bright objects in the five most luminous 
galaxies of the Virgo Cluster Survey. For NGC 3377 and 
NGC 4278 the aperture corrections from a 5 to 10 pixel ra- 
dius were determined with APFILE and found to be —0.132 
in g and —0.170 in z for the former and —0.09 in g and 
—0.13 in z for the latter. For all 14 galaxies the final cor- 



4 This research has made use of the NASA/IPAC 
Extragalactic Database (NED) which is operated by the 
Jet Propulsion Laboratory, California Institute of Technology, 
under contract with the National Aeronautics and Space 
Administration. 



rection from 10 pixel s to infinity was —0.1 in g and —0.12 
m z, from Table 5 of [S irianni et all (|2005h . Extinction cor- 
rections were applied according to A g = 3.634 x E(B — V) 
and A z = 1.485 x E(B-V), where the reddening E{B- V) 
value is from NED and the ratios correspond to the spectral 
en ergy d istribution of a G2 star from T able 14 of|S irianni et 



[al ( |2005| ) and are appropriate for GCs \ Jordan et al.|2004 [ 
The final values A g and A z for each galaxy are listed in 
Table [3j The magnitudes were finally transformed to the 
AB system using as zero poi nts for g and z, 26.068 mag and 
24.862 mag from Table 10 of ISirianni et al N2005I. 



3.4. Comparison with ACSVCS 

We have not used the published catalogues of GCs of 
the ACSVCS (Jordan et al. 20091 as they were not avail- 
able when the present study began. Nevertheless, we com- 
pare the g and z magnitudes of this study with those of 
the ACSVCS for the combined sample of NGC 4486 and 
NGC 4649. The ~ 500 sources of these two galaxies before 
visual inspection were matched with the ACSVCS cata- 
logue. In Figure [3] the difference between the g and z mag- 
nitudes of the present study and those of the ACSVCS are 
shown as function of the g and z magnitudes of the present 
study. The top panels are for the model magnitudes and the 
bottom panels for the av erage correction ap erture magni- 
tudes of the ACSVCS (see | Jordan et al.|2009| ). The median 
difference between the g magnitudes is ~ —0.009 for the 
model and ~ 0.006 for the average correction aperture mag- 
nitudes. On the other hand, the difference between the z 
magnitudes is ~ 0.030 for both comparisons, slightly larger 
for the second case. The scatter is smaller for the comparis- 
son with the average correction aperture magnitudes of the 
ACSVCS, as expected, once the same method of average 
correction obtained from bright GC candidates is adopted 
in both cases. The ACSVCS model magnitudes include the 
correction for the size of each cluster. We conclude that 
either we have underestimated the aperture correction for 
the z band or that the ACSVCS study has overestimated 
it. This difference however, is not expected to affect signif- 
icantly the results of this series of papers. 



3.5. NGC 4365: LIRIS vs. SOFI comparisson 

LIRIS imaging was taken for the galaxy NGC 4365 which 
showed evidence for having intermediate age clusters in pre- 



vious optic al /NIR photometry (with ISAA C /VLT, Puzia 
|et al. 2002[) and was thoroughly studied by Larsen et al 



2005 We~compare d the NTT/SOFI GCs for NGC 4365 



from Larsen et al. (2005) and found 78 objects in com- 
mon, which are shown in Figure |4j Figure [5] shows a com- 
parison of the LIRIS and SOFI K-band photometry, with 
the weighted mean difference being 0.113 ± 0.005 for all 
matched sources. It seems that the difference between the 
LIRIS and SOFI photometry is magnitude dependent as 
the differences become larger as on goes to fainter mag- 
nitudes, with the exception of the brightest magnitudes 
(~17, 18 mag). It is difficult to understand the origin of 
this trend which is oppo site to what is seen for the SOFI 
and ISSAC comparison ( |Larsen et al.j 2005 ) where magni- 
tude differences become smaller with fainter magnitudes. If 
one removes the object with K(LIRIS) ~ 17, the brightest 
of the 5 objects located close to the right border of Fig. 
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Fig. 3. Comparison of the ACS magnitudes of the present 
study and those from the ACSVCS for the combined sample 
of NGC 4486 and NGC 4649. The solid red lines indicate the 
median differences. The top panels are for the model mag- 
nitudes and the bottom panels for the average correction 
aperture magnitudes of the ACSVCS. 



[4j the trend does not appear as significant. The mean and 
median differences between LIRIS and SOFI photometry 
are 0.139 ± 0.031 and 0.094 respectively. When considering 
only sources brighter than 20 mag, the mean and median 
values are 0.09 ± 0.02 and 0.075, respectively, while the 
weighted mean is 0.112 ± 0.005. The di fference found by 
Larsen et aL] ( |2005[ ) between the ISAAC ( |Puzia et al.|2002[ ) 



and SOFI data was 0.102 ± 0.026. Thus, the LIRIS mag- 
nitudes are similar to those from ISAAC, while the SOFI 
measurements show an offset of ~ 0.1 mag. In terms of ages 
and metallicities, this magnitude difference of ~ 0.1 mag in 
the fC-band means that the LIRIS and the ISAAC photo- 
metric systems will yield more metal-rich and younger GCs 
relatively to the SOFI system. This offset, although large, 
is reasonable given the difficulties of accurately calibrating 
NIR photometry due to the bright sky background, uncer- 
tain photometric zero-points, and aperture and photomet- 
ric corrections. For a more detailed discussion see La rsen etl 



al. (2005) 



4. K-band photometric accuracy 

The errors returned by PHOT might not be sufficiently ac- 
curate. They are determined from an image that had the 
galaxy light removed from, and when its individual expo- 
sures were combined, there was subpixel shifting. The com- 
bination of these two effects can be responsible for smooth- 
ing out the original noise fluctuations. When the photo- 
metric errors are estimated they will not account for these 
noise fluctuations or the diffuse galaxy light over the clus- 
ters. Therefore, the output errors could be underestimated. 
To estimate the true uncertainty in the photometry, a cou- 
ple of tests were performed. 
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Fig. 4. WHT/LIRIS K-band image of NGC 4365. The 78 
object s in common with NTT/SPFI from Larsen et al. 
(2005) are marked. North is up and east is left; the grey 
scale is inverted. The cross-shaped feature is a residual from 
the galaxy model subtraction. 
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Fig. 5. Comparison of WHT/LIRIS and NTT/SOFI GC 
magnitudes. The weighted mean difference is 0.113± 0.005, 
and the solid line indicates this value. There were 78 
matches. 



4.1. Artificial star tests 



Artificial star tests with BAOLAB.MKSYNTH (Larsen 



1999) were performed. At these distances, GCs are point 
like sources when observed from the ground. This task gen- 
erates a synthetic image consisting of point sources with a 
shape based on the PSF of the image. The PSF was cre- 
ated with the IRAF PSF and SEEPSF tasks from bright 
stars in the NGC 4278 K band image. A grid of evenly 
spaced sources with the same magnitude was generated for 
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16 < K < 22, binned in 0.5 magnitudes, and added to real 
images of one of the galaxies. This corresponds to a total 
of 110 sources added per magnitude interval, which adds 
up to 1430 sources. The grid was completely regular, and 
the separation between the added sources was 80 pixels. 
The galaxy image used for this was NGC4278 because its 
field is not particularly crowded and there were sufficient 
star-like objects to create a decent PSF. The seeing of this 
galaxy is ~ 1.2", while some other galaxies have better see- 
ing (eg. - 0.9 for NGC4365, NGC4406 and NGC4649). 
With worse seeing the S/N drops, and therefore that the fi- 
nal photometric scatter that we derive will be overestimated 
for the cases with slightly better seeing. Having images with 
point-like sources, we recovered the objects when possible 
and tested the difference between the input and the recov- 
ered magnitudes. The photometry was performed in the 
same way as for the real GC candidates in the galaxies. We 
allow a maximum shift of 1 pixel from the input position 
when recovering the artificial objects and apply a cut in 
the error at K Qn - = 0.5. The last two cuts mean that we 
do not keep objects whose positions shift more than 1 pixel 
from the position where the artificial object was added. In 
Figure |6]the relation K crT , as measured by PHOT, as a fun- 
cion of K is plotted as filled stars for the final GC sample of 
NGC4278, and the output K err as a function of the K out- 
put magnitude for the artificial star test objects is plotted 
with filled points for the criteria: K m < 0.5, x s h < —1, 1 
and y sh < -1, 1. By comparing the GCs of NGC4278 in 
Fig. [6] with the artificial star objects, it is readily visible 
that the relation between K vs. K orl for these two sets of 
objects resemble each other. Thus, applying these 2 criteria 
to the artificial star objects, we are able to reproduce fairly 
well the K and K m curve for NGC 4278. In the upper panel 
of Figure [7] the input K binned in 0.5 magnitude vs. the 
difference between the input K magnitude and the output 
magnitude is shown. The magnitude difference increases as 
we go to fainter magnitudes, and we notice a positive offset 
in the difference around K = 20.5. This offset is attributed 
to recovering a brighter object than it actually is. This can 
happen because when there is a brighter object nearby than 
the input artificial object, PHOT will tend to measure this 
brighter object, even when the maximum shift allowed is 
1 pixel. By applying the x s h < —1,1, y s h < — 1,1 crite- 
ria, these objects with larger position offsets were removed 
and thus the K- m — K out went from the maximum value of 
5.2 to 1.4 for K = 21. Another reason for this to happen 
is that the random fluctuations in the background are so 
high that, for faint objects, a positive fluctuation will make 
the object's estimated magnitude brighter than the input 
magnitude. Also, large negative fluctuations would result 
in a negative flux within the aperture and thus no magni- 
tude would be measured. This is a reason why some ob- 
jects are not detected. A small, negative fluctuation within 
the aperture could also result in a underestimation of the 
magnitude with K in — K out < 0. In the bottom panel, 
the output binned K PHOT error estimate for NGC 4278 
(JCerrN4278) vs. the estimated sigma for the corresponding 
binned magnitude is shown. The best-fit line between these 
two quantities is plotted as a solid line, and the relation is 
given by: 



a = 2.33 x K e 



0.03 



(4) 
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Fig. 6. K err as measured by PHOT as a function of K for 
GCs of NGC 4278 marked as stars and for the artificial star 
test output objects marked with filled dots. 



4.2. SExtractor uncertainty determination 

The accuracy of t he K-band photometry w as also tested 



with SExtractor ( Bertin fc Arnouts 1996[) for one case 



Unlike PHOT, SExtractor QBertin & Arnouts||1996p takes 



into account that there was a galaxy underlying the de- 
tected sources when determining the photometry, which 
means that another source of uncertainty will be added to 
the final photometric error budget. While the photometric 
uncertainties estimated by PHOT are based on the back- 
ground annuli, SExtractor allows the use of two images: one 
with galaxy subtraction and another without it. This latter 
image is used to create a weighting image that is then used 
to calculate the background and thus the errors. As a test 
case we chose NGC 4486. This galaxy has ~300 clusters 
in the LIRIS/ACS field of view, and this will give a more 
robust comparison in terms of number of clusters. We ran 
SExtractor for NGC 4486 and the program detected 77% of 
the 301 sources that fullfill the criteria described in Sect. 5. 

The upper panel of Fig.|8]shows the K magnitude deter- 
mined by PHOT for NGC 4486 vs. the difference between 
the magnitudes determined by SExtractor and PHOT, and 
the weighted mean difference is 0.047 ± 0.008. Note that 
around i^pHOT ft 19 there is an excess of sources with 
a -^Extractor - ifpHOT ~ -0.7. These objects have par- 
ticularly high negative background fluctuations within the 
annuli used to remove the sky contribution, which could be 
the cause of this negative feature. The bottom panel shows 
the magnitude error determined by PHOT for NGC 4486 
vs. the magnitude error determined by SExtractor and the 
best second order polynomial fit between these, given by 
(5), is 



-Kerr SExtractor = —22.88 X K eTT PHOT + 5.40 X K t 



err PHOT 



(5) 



Note that the correlations between the uncertainties in 
the photometry determined by PHOT and the ones deter- 
mined from artificial star tests and by SExtractor, given by 
Eqs. (3) and (4), respectively, are very similar. The con- 
clusion from these tests with BAOLAB.MKSYNTH and 
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0.00 0.05 0.10 0.15 0.20 
K 

rx err N4278 

Fig. 7. Upper panel: input K binned in 0.5 magnitude vs. 
the difference between the input K magnitude and the out- 
put magnitude. Bottom panel: the output binned K PHOT 
error estimate for NGC4278 (Kerr at427s) vs. the estimated 
sigma for the corresponding binned magnitude. The best fit 
between these two quantities is plotted as a solid line. 



0.5 ; 




0.00 0.05 0.10 0.15 0.20 0.25 0.30 

^err PHOT 

Fig. 8. Upper panel: K magnitude determined by PHOT 
for NGC4486 vs. the difference between the magnitudes 
determined by SExtractor and PHOT. Bottom panel: The 
magnitude error determined by PHOT for NGC4486 vs. 
the magnitude error determined by SExtractor. The best 
second order polynomial fit between these two quantities is 
plotted as a solid red line. 



SExtractor is that PHOT errors arc underestimated by 
about a factor of 2, and this has to be taken into account 
in any further analysis. 

5. Selection of globular cluster candidates 

The coordinates of the objects in the ACS z frame were 
transformed to the LIPJS frame with the GEOMAP and 
GEOXYTRAN tasks in IRAF. We used as many objects as 
possible to compute the transformation, from a minimum 
of 10 to a maximum of 25. With GEOMAP the transfor- 
mation from the ACS to the LIRIS coordinate system was 
computed with an output RMS of the transformation for all 
galaxies smaller than 0.25 LIRIS pixels. The ACS z-band 
coordinates were then transformed to LIRIS coordinates 
and aperture photometry was measured with PHOT with 
a centroid algorithm allowing a maximum shift of 1 LIRIS 
pixel. The differences between the coordinate systems of 
the g and z images were a small fraction of a pixel (< 0.1 
ACS pixel). Running photometry in the g image with the 
coordinates from the z image was accurate enough. 

A typical GC has a size of ~ 3pc, and at the distance 
of the galaxies of the sample 1 pixels 4 pc. The spatial res- 
olution of the ACS images allows classification of GCs by 
their spatial extent (they are somewhat more extended than 
stars at the ACS resolution) , and it is thus straight-forward 
to eliminate most of the potential contamination from stars 
and/or galaxies. Effective radii (i? e ff) were measured in g 



images using the ISHAPE code described in Larsen (1999). 
The sizes were measured in these images due to the fact that 
the PSF in this band is more centrally concentrated. King 
models with concentration parameter c (c — r t idai/r C ore), 
fixed to 30 were convolved with an empirical PSF derived 
from a 47Tuc g band ACS image taken around the same 
time the ACSVCS started. The PSF was created with stars 
in this 47Tuc image with the PSF and SEEPSF tasks on 
IRAF. To convert the measured sizes from ISHAPE into 



parsecs, galaxy distances from Tonry et al. ( 2001 ) were used 
for all galaxies with the exception of NGC 4570, whose dis- 
tance was chosen to be the mean of t he available method s 
presented in NED ( [Jordan et al.|2005| and [Tully et al.|1988]) . 
All these distances are listed in Table [I] Only objects with 
a photometric measurement in all three bands (g, z and K) 
were considered to be in the sample. 

In the left panel of Figure [9] we show the optical colour 
magnitude diagrams (CMDs) for the GC candidates de- 
tected in all three bands and that had sizes successfully 
measured, ie., where the KING30 profiles were successfully 
convolved with the empirical PSF. Bimodality is readily 
seen for several galaxies (NGC 4486, NGC 4649, NGC 4526 
and NGC 3377). The NGC 4649 CMD looks different than 
the CMDs for the rest of the galaxies due to several faint, 
blue objects with (g — z) < 0.5, which are likely to be con- 
taminants from a companion spiral galaxy that can be seen 
in the upper right corner of the NGC 4649 image in Figure 
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[2 In the right panel the optical/NIR CMDs for the same 
objects as in the left panel are shown. 

To obtain a cleaner sample of star clusters, only objects 
brighter than g=23 and with optical colours in the range 
0.5 < (g — z) < 2, exp ected for clusters in early- type galax- 
ies ( Peng et al.|2006 ), were kept in the sample. A box with 
dashed lines is shown in the left panel of Figure [9] defin- 
ing the regions that fall within these colour and magnitude 
cuts. The objects that fall inside the box in the left panel are 
indicated by a different colour/symbol in the right panel. 
We define a lower magnitude cut in an optical colour due 
to the larger errors in K, specially at fainter magnitudes. 
The magnitude cut translates roughly into an error cut, as 
can be seen in Figure [10] This Figure shows the PHOT 
errors for g, z and K . Objects brighter than g=23 have g 
and z errors of about 0.05 mag while the bulk of objects 
have K-band errors below 0.2 mag. Following the ISHAPE 
measurements, we added a size range criterion. The objects 
had to have sizes between 1 < R c ff(pc) < 15 to be kept in 
the final sample of GCs. This is the range that includes the 
great majority of Milky Way GCs. The lower cut in sizes of 
Resipc) > 1 removes objects that are essentially unresolved 
and foreground stars whereas the upper cut of R c g (pc) = 15 
excludes background galaxies. With these cuts one should 
keep essentially star clusters. Although some very distant 
background galaxies might still be contaminating the sam- 
ple, we did not estimate the contamination from these. To 
illustrate the size criteria, we show the distribution of effec- 
tive radii (R c g) for the NGC4486 GC candidates in Figure 
12 that had already satisfied the colour and magnitude cri- 
teria: 0.5 < (g — z) < 2.0 and g < 23. Note the peak around 
R e s ~ 2.5, which is expected and is in accordance with 



Jordan et al.|(2005|. By adopting 1 < R e g < 15pc as crite- 



rion ~ 8% of the candidates that had sizes successfully mea- 
sured for NGC 4486 drop out of the sample. In Figurc[TT]thc 
g magnitudes are shown as function of R c ff for the GC can- 
didates that were detected in the three bands and that had 
sizes successfully measured, as in Fig.[9j The final sample of 
GCs is indicated by a different colour /symbol. Most GCs in 
the final sample have R e g ~ 3. Note the clear separation be- 
tween point-like sources (stars) and GC candidates, < 1 pc. 
Moreover, we excluded all objects in NGC 4649 that were 
located inside a rectangle of size 40" x 90" that contains 
the smaller companion galaxy. 

Furthermore, since the resolution of ACS is much su- 
perior to that of LIRIS, a final criterion for GC selection 
was added. If GC candidates are too close to each other 
in the ACS images they have a chance of appearing as one 
bigger source in the LIRIS K-band image, i.e they might be 
blended with another source in the K-band. The better the 
seeing, the smaller this chance. To take this into account the 
cluster candidates were visually inspected and the ones that 
looked like they might be blended with another source were 
flagged. We defined a minimum distance between neigh- 
bouring sources not to be blended with each other as -Rbiend- 
We also assigned a criterion for the brightness of the neigh- 
bouring object. The neighbouring source should be at a 
distance < -Rbicnd and the flux of this neighbour should be 
> 10% of the flux of the GC candidate for it to be consid- 
ered a blend. For example, for a GC candidate of NGC 4486 
that was flagged as a possible blend, the separation between 
the neighbour and the GC candidate should be < 30 pix 
in the ACS image (~ 1.5") and the flux of this neighbour 
should be >10% of the flux of the GC candidate for it to be 
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Fig. 10. PHOT errors in g, z and K as a function of g 
magnitude for NGC 4486 as an example. The dashed line 
defines the cut, only objects brighter than g=23 were kept 
in the sample. 



removed from the sample. i?bicnd was determined by visual 
inspection of the smallest distance between two sources in 
ACS to appear as two sources in LIRIS. i?bicnd is related to 
the seeing: if seeing = 0.8—1.0", -Rbicnd = 25 ACS P i X ; if see- 
ing = 1-1.2", i?bic„d = 30 ACS plx ; if seeing = 1.3-1.5", 
^bicnd = 35 ACSpix. -Rbicnd an d the number of blended 
sources (Nbiend) for each galaxy are listed in Table [3j 
Together with eliminating the blends, we checked and ex- 
cluded those GC candidates that were obvious background 
galaxies. Also, objects that fell in the residuals of the K- 
band model subtraction in the very centre (~ 24") of the 
galaxy where removed. The fraction of eliminated objects 
in the visual inspection ranged from 3 — 20%, and this de- 
pends on how cluster rich a galaxy is. For the majority of 
the cases, 10% of the cluster candidates were eliminated; as 
expected, the richest cluster galaxy NGC 4486 had 20% of 
the candidates removed from the sample due to crowding. 
In Fig. [13] we show the final reduced and model subtracted 
NGC 4486 image where we ran photometry on. GCs that 
fullfill the criteria defined above are shown as circles. In 
Table [4] the g, z and K photometry, astrometry and size 
information is given for a sample of these NGC 4486 GCs 



shown in Figure 13 



In summary, the selection criteria are: 



— magnitude cut g < 23. 

— optical colour cut 0.5 < (g — z) 



< 2. 
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Fig. 9. Le/£ panel: Optical ((5 — z) vs. g) colour magnitude diagrams for globular cluster candidates detected in g, z and 
K that had a size assigned to by ISHAPE in grey. The box with dashed red lines defines the colour and magnitude cuts. 
Only objects brighter than g=23 and with optical colours in the range 0.5 < (g — z) < 2 were maintained in the sample. 
The final sample of GCs after visual inspection and size cut (see text for details) are indicated as blue squares. Right 
panel: Optical/NIR ((g — k)vs.g) colour magnitude diagrams for the same objects as in the left panel. Objects inside 
the box with dashed lines that are grey dots in the left panel are shown as red triangles in the right panel. Blue squares 
are as in the left panel. 



— size cut 1 < R c g < 15 eliminating background galaxies 
and foreground stars. 

— magnitude uncertainty cut K crr < 0.5. 

— visual inspection, flagging blends (see text for details) 
and still left over spurious objects such as obvious back- 



ground galaxies and objects over the residuals of the 
K-band model subtraction in the centre of the galaxy. 
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Galaxy 


ID 


RA(J2000) 


DEC(J2000) 


a 




z 


az 


K 


<tK 


-Rcff(pc) 


NGC 4486 


1 


187.71700 


12.36329 


22.944 


0.014 


21.732 


0.012 


19.462 


0.094 


2.994 


NGC 4486 


2 


187.69138 


12.36336 


22.572 


0.010 


21.242 


0.008 


19.369 


0.059 


2.599 


NGC 4486 


3 


187.72313 


12.36383 


22.384 


0.010 


20.935 


0.007 


18.662 


0.025 


2.373 


NGC 4486 


4 


187.69490 


12.36397 


22.560 


0.011 


21.028 


0.007 


18.697 


0.023 


2.316 


NGC 4486 


5 


187.71342 


12.36401 


22.947 


0.013 


21.436 


0.010 


19.095 


0.053 


1.412 



Table 4. ACS and LIRIS photometry and sizes for the GC system of NGC 4486. The full version of this table, containing 
the other 13 GC systems presented in this study is available online at the CDS or upon request from the authors. 
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Fig. 11. The g magnitude as a function of R c ff . The box 
with dashed lines defines the size cut: 1 < R c g(pc) < 15 and 
the magnitude cut g < 23. The final sample, after visual 
inspection is indicated as blue squares, as in Fig. [9| 
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Fig. 12. The distribution of sizes for the GC candidates of 
NGC 4486 as an ilustration after the colour and magnitude 
criteria were applied: 0.5 < (g — z) < 2.0 and g < 23. 



6. Summary 

In this paper, we have introduced a homogeneous survey 
of globular cluster systems in early-type galaxies from a 
combination of NIR and optical photometry. We have ac- 
quired deep K-band images for 14 early-type galaxies with 
the LIRIS instrument on the William Herschel Telescope 
and combined with the HST/ACS optical g (F475W) and z 
(F850LP) bands. Based on sizes, colours and optical bright- 
ness we have defined a sample of GCs for each galaxy. We 
have described the observing strategy and the data reduc- 
tion procedures. Compar ing the LIRIS data to NTT/SOFI 
data ( |Larsen et aL|2005| > for the galaxy NGC 4365 we find 
a difference of ~ 0.1 mag showing that our photometry is 
reliab l e and agrees better with the ISAAC results by |Puzia] 



et al.| ( |2002[ ). The errors from PHOT for the GCs are found 
to be underestimated by roughly a factor of two. 

This survey has the great advantage of being homoge- 
neous: data from the same instrument, with roughly the 
same exposure time, reduced and analysed in the same 
way. It is also aided by homogeneous optical ACS/HST 
data, which helps the decontamination of the GC sample 
by its resolution, which allows size measurements. In total, 
we have assembled a sample of ~ 1300 GCs in 14 galax- 
ies. For the majority of the galaxies we detect ~ 70 GCs 
with g < 23. NGC 4486 and NGC 4649, the cluster-richest 
galaxies in the sample contain 301 and 167 GCs, respec- 
tively. These are the largest number of GCs so far for which 
NIR photometry has been obtained. 
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Fig. 13. Final reduced and model-subtracted NGC4486 
image with its final optical/NIR GCs as described above. 
North is up and east is left; the image scale is inverted. 



In the upcoming papers of the series we will address the 
science results of the data presented here. 
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